Falsification of the MRI in the Two-Dimensional Axial Symmetry 
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We analyze a linear perturbation scheme for a two-dimensional background plasma, which is 
rotating with a differential frequency and is embedded in a poloidal magnetic field. The plasma 
profile is axially symmetric, both in the background and the perturbation approximation and the 
azimuthal magnetic field is requested to vanish identically. The main assumption of the present 
study is that in the non-stationary regime still holds the dependence of the angular frequency on 
the magnetic surface function only, which, on the steady background equilibrium, is ensured by the 
validity of the co-rotation theorem (Ferraro 1937). Indeed, such a restriction of the model is rather 
natural and it implies that the azimuthal component of the linear plasma shift is reabsorbed in the 
expression for the non stationary electric field (in principle, at any order of approximation) and can 
no longer provide a non-zero azimuthal component of the magnetic tension field. As a result, the 
Magneto-rotational instability is suppressed and the magnetic field has the effect to stabilize the 
plasma configuration with respect to the pure hydrodynamical case. 

Keywords: Plasma physics; Axial symmetry; Magneto-rotational instability (MRI); Co-rotation theorem 



The presence of rotation in the quasi-stationary state 
of a magnetized plasma is a common feature, both in as- 
trophysical and laboratory systems. In stellar accretion 
disks, the plasma rotates because of its orbital motion in 
the gravitational field of the central body [l|. In Toka- 
mak experiments, the emergence of rotation profiles is 
observed as spontaneous tendency of the system evolu- 
tion and its origin is still questioned, though three main 
proposals stand Q. 

The coupling between the magnetic field, in which 
the plasma is embedded, and the differential rotation 
of its different layers, lead to the appearance of unsta- 
ble modes, named Magneto-rotational instability (MRI), 
first discovered by Vclikhov in 1959 Q and applied to an 
astrophysical context by Chandrasekhar in 1960 Q. This 
effect is relevant in stellar accretion disks, as well as in 
Tokamak configurations in order to trigger the onset of a 
turbulent behavior in both these scenarios of plasma dy- 
namics. The merit of having properly characterized this 
crucial role of MRI in determining the transition from 
a laminar to a disordered regime of a differentially ro- 
tating plasma is due to Balbus and collaborators in the 
early nineteens @ (see also @). 

In the case of a plasma disk, surrounding a compact 
astrophysical object, the presence of the MRI is the only 
viable mechanism to generate a turbulent transport, re- 
sponsible for accretion, when restated as a laminar vis- 
corcsistive behavior, according to the original idea for 
angular momentum transport [7|. 

In the Tokamak plasma profiles are commonly regis- 
tered collective behaviors of the configuration, associated 
to both toroidal and poloidal rotation velocities. Such 
a phenomenon, dubbed spontaneous rotation, is often 
favourable to the plasma confinement, allowing steepen 
profiles of temperature and density toward the torus axis. 
However the plasma rotation is naturally suppressed in 



favor of a turbulent evolution and, often, it takes place a 
real interplay between the laminar and chaotic regimes. 
In such a interchange, the MRI is expected to play a role 
in originating the unstable modes, then evolving toward 
the onset of turbulence, whose associated transport is, in 
gneral, a source of spread for the plasma confinement 

For a discussion on the parallelism existing between 
the axial symmetry of a stellar accretion disk and a Toka- 
mak configuration see @ , which was recently enforced by 
the development of local disk profiles [HI El character- 
ized by magnetic field micro-structures resembling well- 
known features of the laboratory plasma physics. For the 
generalization of such profiles to the global radial equi- 
librium of a steady state of an accretion disk see |13 |. 
while their impact on the implementation of the Shakura 
prescription in a fully magneto-hydrodynamical scheme 
is analyzed in [l3j . 

Here, we consider a two-dimensional axisymmctric 
configuration of the plasma, built up by a steady back- 
ground, characterized by a rotating plasma, and a first 
order linear perturbation theory, which introduces a time 
dependence in the model. The normal modes analysis 
is performed around a local configuration of the back- 
ground, determined by fiducial values of the radial and 
vertical coordinates. The analyzed configuration is re- 
stricted by requiring that the angular velocity is, at any 
order of approximation, a function of the magnetic sur- 
face only. For the stationary axisymmetric configuration 
this feature is assured by the co-rotation (also known as 
iso- rotation) theorem of Ferraro (HI, but, in the linear 
time dependent perturbation scheme, it plays the role 
of a specific assumption on the considered configuration. 
Furthermore, we require that the azimuthal component 
of the magnetic field identically vanishes. Its exact va- 
lidity is crucial in the derivation below, however such a 
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restriction could be motivated a posteriori both for as- 
trophysical and Tokamak systems. 

Starting from these two statements, we arrive to 
demonstrate that, in such a system, the MRI is sup- 
pressed and that a sufficiently strong magnetic field, is 
even able to remove pure hydrodynamical instabilities. 

The mechanism responsible for the MRI suppression is 
the possibility of reabsorbing the azimuthal component 
of the plasma shift into the poloidal electric field. Indeed, 
we are in practice working in the Coulomb gauge, deal- 
ing with a time-dependent electric potential, made up of 
the magnetic surface function only. Thus, the time de- 
pendence of the perturbed magnetic field arising in the 
plasma is fixed by the pure poloidal component of the 
shift and the azimuthal component of the magnetic ten- 
sion vanishes, implying the MRI suppression. 

This result is of impact both in the astrophysi- 
cal and laboratory settings, because it shows how, if 
the azimuthal magnetic field is sufficiently small and 
the angular frequency is forced to strictly follow the 
magnetic surface structure, the fundamental Magncto- 
hydrodynamical instability, affecting the plasma rota- 
tion profiles, i.e. the MRI paradigm is naturally re- 
moved. Indeed in the Tokamak configurations the tran- 
sition from spontaneous rotation to turbulent regimes is 
expected to take place by non-linear transitions in re- 
gions of the system parameters (background magnetic 
field and perturbation wavelength), where MRI is any- 
way suppressed 0, [H, [f5l - [l8j . Otherwise the issue of the 
present analysis strongly concerns the angular momen- 
tum transport in astrophysics (see the discussion at the 
end of the paper). In fact, when the MRI is absent in the 
linear backreaction of an accreting disk plasma, there ex- 
ist no firm alternative scenario to generate the effective 
dissipation required by the standard model for angular 
momentum balance across the configuration [l|, 0, [r| ■ 

Having in mind the small scale linear perturbations to 
a purely rotating plasma disk, let us consider the system 
made of the Faraday law and of the electron force balance 



d t B = -cV A E, E+-AB = 0, 
c 



(1) 



where E and B denote the electric and the magnetic field 
respectively. Equation (p} admits the solution 



E = -V$ - -d t A. 

c 



(2) 



where A is the vector potential, such that the magnetic 
field B — VAA, while $ denotes the time dependent elec- 
tric potential (indeed, we are operating in the Coulomb 
gauge V • A = 0). Substituting this solution into the 
electron force balance, we get the fundamental relation 



V$+ -d t A = - A B 

c c 



(3) 



We now specialize this scheme to the two-dimensional 
axial symmetry, by assigning a vector potential of the 



form A = (ip(t, r, z 2 )/r)e<p, where we adopted cylindrical 
coordinates {r, </>, z} and we take the magnetic flux sur- 
face function ifj as depending on z 2 because the reflection 
symmetry with respect to the equatorial plane z = is 
postulated (for a stellar accretion disk, this is due to the 
symmetry of the gravostatic equilibrium) . The resulting 
magnetic field takes the explicit form 



B = — d z 4>e r H — d r ipe z 
r r 



(4) 



We also consider the following velocity field in the plasma 



v = w(-0)re0 + v p , 



(5) 



where v p = v r e r + v z e z denotes the poloidal velocity field 
in the meridian plane. This choice of the magnetic and 
velocity fields contains two main assumptions, the van- 
ishing nature of the toroidal component B^ and that the 
angular frequency of the disk depends, on the magnetic 
flux surface, u> = ui{ip). Nonetheless, the considered sys- 
tem is rather general in axial symmetry and it does not 
contain, in principle, any restriction able to suppress the 
MRI. Using the expressions above for B and v to rewrite 
the relation ([3]), we get 



V$ + — = — -^-V^ [v p ■ Vip 

c r c cr x 



(6) 



From the poloidal and azimuthal components of this 
equation, we arrive to write down the two key relations 



dip c 
d t ip = —v p ■ V%1>. 

In what follows, it will be convenient to use relation 
in its vector form 



d t A 



AB. 



(7) 

(8) 
© 

(9) 



We now consider a linear perturbation scheme to the 
background equilibrium of a rotating plasma disk, hav- 
ing angular velocity lo = cjo(V'o( r o, z o))j a mass density 
p = po(ro,ZQ) and the pressure p = Po( r cb ^o)' an< ^ t> e " 
ing embedded in a magnetic field described by a surface 
function V'o^Oj z o) (f° r a detailed discussion of the static 
force balance, see the discussion below, at the end of the 
linear perturbation calculus). 

Here ro and zq are two fiducial values of the radial 
and vertical coordinates respectively, around which we 
develop the local behavior of the perturbations, whose 
wavelength is assumed to be sufficiently small to ex- 
plore only the region nearby such background circumfer- 
ence. In this approximation limit, we take the magnetic 
surfaces in the form -0 = iPo{i~o>Zq) + ^i(t,r,z), with 
I ipi ipo I and 



(V^i) |<| (VoV'o) 



(10) 



where Vo denotes the gradient operator with respect 
to the fixed space coordinates ro,ZQ. Such a condition 
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clearly states that the perturbed magnetic field is much 
smaller than the background one, as requested by the lin- 
ear approximation. However, we also require the natural 
inverse relation for the second order gradients (current 
density), i.e. 

| (V 2 fa) »| (V^o) I • (11) 

This paradigm for the gradient fields corresponds to the 
natural hierarchy relation in the plasma. We introduce 
the poloidal shift vector £ p of the plasma, by defin- 
ing the corresponding poloidal velocity as v p (t, r, z) = 
dt£ p (t, r, z). While, the angular velocity admits the local 
development 

u = u(fa(r Q - £ P r, z - £, pz ) + fa) ~ 

w - u'orZpr ~ u' 0z £ pz + LO fa , (12) 

where we introduced the notation 

' duj\ ( g?u/ 



J 0x 



dx 



,w 



(13) 



x=x \—r / i/)=tj> 

From this expression, it is immediate to reach the iden- 
tification dt£$ = wofa- In order to focus our analysis to 
the main features of the MRI, we consider incompress- 
ible perturbations, i.e. we require V • V\ = and hence 
we get V • £ p = 0. By equation (0) we get the following 
expression for the perturbed magnetic field 



Bi = (Bo ■ V)£ P 



(14) 



Thus, the magnetic tension force f v m ) takes the form 
(the background pressure and tension are negligible by 
means of the inequality (fTTj) ) 



f< m > = p v A (B -V)%, 



(15) 



here Bo is the unit vector in the magnetic field direc- 
tion and v A denotes the Alfven speed associated to the 
background. Taking the Fourier representation of the 
perturbation £ p 



Cp = Zp(i) exp{i(k p ■ f p } . 



(16) 



where the index p denotes the poloidal component of the 
wave and position vectors. We note that using condi- 
tions (fT0|) and (fTTj) we take k p r 3> 1 and k p r <C 1^0/ fa- 
in order to select pure Alfven perturbations, i.e. to elimi- 
nate the magnetic pressure from the problem, we restrict 
our attention to wave- vectors of the form k p = kBo, so 
that, as a consequence of the incompressibility (namely 
k p ■ £ p = 0), we obtain the constraint Bo ■ £ p = 0. Then, 
the magnetic tension takes the form 



f M = -pok 2 v% = -po^Up . 



(17) 



The momentum conservation equations, for the velocity 
perturbations Vi — V p i + i^ie^, take the form 

po (9 t u r i-2cj ^i)=T r (m) ,(18) 

Pod t v zl =Ti m \ (19) 

po (dtv^i + 2oj v r i + rQUj' 0r v r i+roU) 0z Vzi) = T^ m) (20) 



e.g. [l|. By means of the natural relations 

v r i = d t £,pr , v z i = dt£ pz , (21) 
v</,i = r (uofa - u' 0r £pr " u' Qz £ pz ) (22) 

the explicit form of the magnetic tension, the momentum 
conservation equations rewrite as 

dt£pr - 2u r uofa =-{Vr + u A )£, P r - Vztpz, (23) 

d 2 t tp Z = -u? A tp Z , (24) 
r w <9 t V'i + 2u d t £pr = , (25) 

where we adopted the notation y p = 2roUJo^'op- H ere i we 
retained the full dependence in £ p and then the partial 
time derivative, since we did not yet Fourier expand fa. 
The dynamics of the vertical shift £ pz is fixed by equation 
(|2~4")) and it yields a stable behavior of the plasma along 
such a direction. Equation (|25"|) can be easily integrated, 
getting the relation 



rouofa 



(26) 



which substituted into equation (|23| gives the following 
equation for £, pr 

dt£pr = -(K-o + UJ A )ipr-yzlpz , (27) 



where /Co denotes the epicyclic frequency, i.e. 
1 

3 
' o 



(28) 



Deriving equation (|27[) two times in t and using equa- 
tion (l23l) we obtain: 



d% r = -(/Cg + ufo%lrr+V*? A ip* ■ (29) 
Using again equation (|27p in equation (j2"9")l one has 
d% r + [Kl + 2cu 2 A )d 2 t i pr +u: 2 A {JCl + u 2 A % r = . (30) 

Taking the following plane-wave representation for the 
time dependence of £ pr 

£pr = ipr exp{-i^i}, (31) 

equation (j3T))) reduces to 

O 4 - (fC 2 Q +2u 2 A )n 2 + u 2 A (IC 2 Q + uj 2 a ) = 0, (32) 

admitting the following solutions 

fl 2 =u 2 A and n 2 = JCl + uj 2 a . (33) 

In the absence of the magnetic field, i.e. u) A = 0, the 
dispersion relations (j3"3")l reduce to the well-known sta- 
bility condition fC 2 > (see @ for a discussion on the 
Rayleigh criterion for local stability). Such a condition 
states that the specific angular momentum of the plasma 
radially increases and it holds as far as the angular fre- 
quency sufficiently smooth decays. When the plasma 
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is magnetized, for having stability, the relation (|33|) re- 
quires K-q + uj 2 a > . Since w\ is intrinsically positive by 
its definition, the effect of the magnetic field consists of 
the stabilization of an unstable part of the spectrum for 
the non-magnetized (neutral fluid) configuration. In this 
scheme, the MRI is completely suppressed and replaced 
by an opposite effect, which stabilizes the differential ro- 
tation profile of the plasma and prevents the onset of a 
turbulent behavior. 

Let us now analyze if the result above remains valid 
in the case of a non-vanishing azimuthal component of 
the background magnetic field, i.e. Bq^ ^ 0. Since the 
(perturbed) poloidal velocity is a divergenceless field, it 
can be expressed via a scalar function 0(r, z) as follows 

ve _ 

v p = — A e . (34) 

Then, in the right-hand side of equation (|6]) we would 
have the additional term S/QB^/r and the possibility 
to reabsorb this contribution into the electric potential, 
requires that it be a pure gradient. Thus, we get the fol- 
lowing constraint for the azimuthal magnetic field compo- 
nent B,f, = rf(Q) , / denoting a generic functional form. 
To understand if such a form of B$ is allowed by the 
considered plasma configuration, we have to analyze the 
background profile. The magneto-static configuration is 
fixed by the basic equation 

V po - poOJ%r e r = F QL - p V $ g , (35) 

where Ft, — (V A B) A B/4tt denotes the Lorentz force, 
$ g is the gravitational potential (it concerns the stellar 
accretion disk configuration only) and po is the back- 
ground pressure. By the axial symmetry restriction, the 
consistency of the azimuthal component of equation (|35[) 



clearly implies Fql,/, = 0. As it is well-known and it is eas- 
ily checked by flD, this condition implies B 0r p = I(tp )/ro, 
here J is a generic function of the magnetic surfaces. As- 
suming this form as valid at any order, a non-zero az- 
imuthal component of the magnetic field is allowed in 
the present scheme, as far as the condition = g(I/r 2 ) 
(with g = / _1 and g{I(ipo)/ r o) ~ 0) holds. In the local 
model we arc considering, any function g having at least 
one zero is suitable for expressing the first order quantity 
0. It is worth noting that, since ip is here taken even 
in z, the radial velocity is odd and vanishes on the plane 
z = Zq, so that locally we have no real accretion toward 
the axis or ejection outward of it. Furthermore even if 
we consider an azimuthal perturbation of the magnetic 
field as far as its background value is zero, such pertur- 
bations are frozen (see eq. ([5])). Anyway, the possibility 
to neglect the azimuthal magnetic field component can 
also relay on phenomenological considerations. In Toka- 
mak configuration this component is much smaller than 
the toroidal (here the vertical) one, almost of the same 
amount of the ratio between the small and large radius 
of the torus, although its presence has a physical role in 
the plasma confinement profile. In the stellar accretion 
disks, the azimuthal magnetic field is essentially due to 
the dynamo effect. Thus, it is associated to the turbu- 
lent behavior of the plasma, whose onset is triggered by 
MRI [2(j. Thus, in both these cases we can regard the 
azimuthal field, if present with a generic form, as a weak 
modification, affecting the picture depicted above. 
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